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The versatile technology of cold atoms confined in optical lattices allows the creation of a vast 
number of lattice geometries and interactions, providing a promising platform for emulating var- 
ious lattice models. This opens the possibility of letting nature take care of sign problems and 
real time evolution in carefully prepared situations. Up to now, experimentalists have succeeded 
to implement several types of Hubbard models considered by condensed matter theorists. In this 
proceeding, we discuss the possibility of extending this effort to lattice gauge theory. We report 
recent efforts to establish the strong coupling equivalence between the Fermi Hubbard model and 
SU(2) pure gauge theory in 2+1 dimensions by standard determinantal methods developed by 
Robert Sugar and collaborators. We discuss the possibility of using dipolar molecules and ex- 
\Q ternal fields to build models where the equivalence holds beyond the leading order in the strong 

coupling expansion. 
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1. Introduction 

New experimental techniques to confine cold atoms in optical lattices have allowed the creation 
of many lattice geometries and interactions. Trapping polarizable atoms or molecules in a periodic 
potential created by crossed counterpropagating laser beams has been an area of intense activity in 
recent years. It is now possible to physically build lattice systems where the average number of 
particles per site and their tunneling between neighbor sites can be adjusted experimentally. This 
has opened the possibility of engineering experimental setups that mimic lattice Hamiltonians used 
by theorists and to follow their real time evolution. Up to now, experimentalists have succeeded 
to implement several types of Hubbard models. Can these successes be extended to lattice gauge 
theory? 

There are problems that are very difficult to address by conventional (digital, classical) compu- 
tational methods but could be dealt with if lattice gauge models could be studied from approximate 
implementations on optical lattices: for instance the phase diagram of QCD at finite temperature 
and significant chemical potential, real time evolution and string breaking. Beyond its success in 
particle physics, quantum gauge theory also have played an important role in the development of 
quantum computing. Decoherence can be an important obstacle for quantum computation. One of 
the first applications of Shor's idea of fault-tolerant quantum computation [1] was Kitaev's toric 
code [2] where the two stabilizers (the plaquette and the cross) are operators that were introduced 
two decades before in the Hamiltonian formulation of lattice gauge theory [3]. For these reasons, 
it seems that some effort should be put in finding optical lattice realizations of quantum gauge 
theories. 

It has been suggested [4] that the approximate emergent local SU(2) symmetry [5, 6] of the 
Hubbard model at strong coupling could be a good starting point for this enterprise. This is briefly 
reviewed in Sec. 2. However, the correspondence Hubbard ~ Heisenberg ~ Hamiltonian Lattice 
Gauge Theory only works up to second order in the strong coupling expansion and there is no 
obvious plaquette interactions. The fermions coupled to the gauge fields have color indices but no 
Dirac spin indices. The strong coupling expansion could be used as a guide to add interactions or 
combine different types of atoms or molecules in order to build Hubbard-like model corresponding 
more accurately to existing lattice gauge models. Some elements of the expanding toolbox that can 
be used for this purpose are provided in Sec. 3. Finally, one should use computational methods 
developed in lattice gauge theory to calculate determinants that have been adapted for Hubbard 
models [7] to test the validity of the approach. 

2. Approximate local symmetry, mesons and baryons in Hubbard models 

The Hubbard model Hamiltonian is 

N 

(i,j},a i'=l 

where t characterizes the tunneling between nearest neighboor sites and U controls the onsite 
Coulomb repulsion. 
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Figure 1: Graphical representation of the interations of the Hubbard model. 

In the strong coupling limit (U 3> t) and at half-filling, which has one fermion per site, the 
previous model is equivalent to Heisenberg model 



<'j> 

with / = 4t 2 /U. The spin operator Si can be written as 
The corresponding Heisenberg Hamiltonian becomes 

<ij> <ij> - 

A constraint must be imposed in order to recover the original Heisenberg model: fLfia = 1- 



(2.2) 



(2.3) 



(2.4) 
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Figure 2: Graphical representation of the half-filling, strong onsite repulsion situation. 
After a particle -hole transformation in the spin down operator [5, 6, 8], 

Jv.fi fuJl < 25 ^,2 

The Heisenberg Hamiltonian can be written exactly as follows 



(2.5) 



(2.6) 
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where the "meson" and "baryon" operators are 

M X = £ % a ^x,a (2.7) 



a=l,2 

and 

/?,.= V 

2 



S.v= I ^^.vA fl =^,,l^.v,2 (2-8) 
0=1,2 

The anti-ferromagnetic Heisenberg model can also be obtained from the strong coupling limit 
of the Kogut-Susskind Hamiltonian of SU (2) lattice gauge theory [6] which is denned as 



H = ~ £ £5£fi + *- £ (V'iAIV'x+M, - *.c.) (2.9) 

with 7= 16/3g 2 °c j3. 



3/ X1 = 2 

«,x,i a, o.x.i 



The approximate gauge invariance allows to show the equivalence between two apparently 
very different mean field theory of the Hubbard models (d-wave super conductivity and "flux 
phase") [5]. The above approximate equivalence has been proposed [4] as a starting point to imple- 
ment lattice gauge theory on optical lattice. However, the equivalence with H in Eq. (2.9) is only 
valid under second order perturbation theory. There are no plaquette interactions for the above 
lattice gauge Hamiltonian and there are no spin indices for the gauge fields. Staggered fermion 
Hamiltonian can be realized by doubling the lattice spacing [8]. Chemical potential generally in- 
troduces sign problem [9] which makes physical simulations desirable. 

3. A toolbox to engineer plaquette interactions 

In the following, we discuss some tools that can be used to build new interactions in optical 
lattice systems. 

3.1 Dipolar fermions in external fields 

External fields can be used to get close to Feshbach resonances or induce dipole interactions 
across sites. From a theoretical point of view, new couplings and long-range interactions can be 
varied to obtain rich phase diagrams [10]. The orientation of the dipole moments can also be ma- 
nipulated to tune the dipole interaction. Considering single species fermions on a two-dimensional 
half-filled square lattice in the xy-plane and fixing the azimuthal angle 0=0, the variation of 
the polar angle 6 has been found [11] to lead to three different phases, as shown in Figure 3. 
For 6 < 6i , the repulsive nearest- neighbor interaction leads to a conventional density-wave phase, 
here denoted as CDW (charge-density-wave) in analogy to condensed matter systems where the 
fermions are charged. With increasing 6 (6i < 6 < 62), nearest-neighbor interactions along the x- 
direction become attractive and a density-wave with high-probabilities along the links sets in. This 
phase has been referred to as a bond-order solid (BOS), or a p-wave CDW. For 6 > 62, the ground 
state has ;?-wave BCS pairing. Generally, with long-range interactions, ordered states with non- 
zero angular momentum becomes favored, possibly involving bond orders. For spin-1 /2 fermions, 
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a p-wave spin-density wave phase has been found for dipolar fermions on a half-filled square lattice 
[12], corresponding to a vector order parameter on the bonds. 
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Figure 3: Phase diagram for spinless dipolar fermions on a square lattice at half-filling as a function of the 
polar angle 6 for the direction of the dipole moments (</> = 0). The CDW phase corresponds to checkerboard 
order, the BOW has different probabilities on the links, and the p-wave BCS has the symmetry of the order 
parameter as shown on the upper-right corner. 



3.2 Combination of different types of atoms or molecules 

An interesting possibility [14] to generate new interactions is to create a bi-partite lattice with 
S-orbits and P-orbits arranged on a checkerboard pattern and tune the relative heights in order to 
have the different orbitals at approximately the same level. This situation can be used to generate 
ring currents which share some feature of the plaquette interactions [15, 16]. Another possibility 
is to use two lattices: one lattice having atoms or molecules that can hop and induce the new 
interactions on the other lattice. Cheng Chin's group uses 6 Li atoms coupled to l33 Cs. Using a 
Fermi-Bose mixture of ultra-cold atoms in a optical lattice, one can couple bosonic gauge fields 
to fermionic matter [18]. Models have been constructed by combining quantum link models with 
spinless staggered fermions [17]. U(l) lattice gauge models could be implemented by using one 
fermion species and two boson species which represent gauge fields. 

3.3 Effective plaquette interactions 

It is clear that higher order terms in the strong coupling expansion will generate plaquette in- 
teractions. More generally, these interactions should appear generically after blocking the Hamil- 
tonian of Eq. 2.6. Recently, an old preprint [20] explaining how non-Abelian gauge interactions 
originate as effective dynamics in models of hopping particles has been made available on arXiv. 

4. Work in progress 

We have been working in several directions. We started up with the 2 dimensional Hubbard 
model and 2+1 dimensional lattice gauge model with 2 colors. Determinantal calculations to 
check the approximate equivalence between the two models are in progress. We are writing a 
Hubbard model code by using the Hybrid Monte Carlo (HMC) algorithm. We plan to add different 
boson species to the existing code in order to compare with models of Fermion-Boson mixtures 
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confined in the optical lattices. Progress also needs to be made in higher orders of the strong 
coupling expansion from Eq. 2.6 to see how effective plaquette terms can emerge. We also need to 
understand the correspondence between the physical temperature used in real experiments and the 
coupling (j8 or g) used in the lattice gauge model with specific 2+ 1 geometries. 

This work was stimulated by the workshop "Critical Behavior of Lattice Models" in summer 
2010 at the Aspen Center for Physics. This work was funded in part by the Mathematical and Phys- 
ical Sciences Funding Program of the University of Iowa. SWT acknowledges support from NSF 
under grant DMR- 0847801 and from the UC-Lab FRP under award number 09-LR-05-1 18602. 
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